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Abstract — Non-minimum phase characteristic of hydro turbine 
shows an opposite initial power surge in the event of frequency 
disturbance, thus possesses widely different characteristic than 
thermal generating unit. This paper presents Automatic 
Generation Control (AGC) of an interconnected two-area 
multiple-units hydro-hydro power system in restructured 
electricity market. The step load perturbation to such a system 
results in heavy frequency oscillations and the system is 
unable to regain its stable state as the positive real parts of 
some eigenvalue pairs confirm the inherent dynamically 
unstable characteristic of the system. To stabilize the 
frequency oscillations, impacts of Superconducting Magnetic 
Energy Storage (SMES) placed at terminal of area and 
Thyristor Controlled Phase Shifter (TCPS) located in series 
with tie-line have been investigated. Two cases such as (a) 
SMES-SMES coordination and (b) TCPS-SMES coordination 
are evaluated to compare their effectiveness to suppress the 
frequency oscillations. The parameters of TCPS, SMES and 
integral gains of AGC loop are optimized through craziness- 
based particle swarm optimization algorithm in order to have 
the optimal transient response of the system under different 
PoolCo and bilateral transaction in restructured electricity 
market. 

Index Terms— AGC, Deregulation, SMES, TCPS 

I. Introduction 

Automatic generation control has been used for several 
years to meet the objective of maintaining the system 
frequency at nominal value and the net tie line power 
interchange from different areas at their scheduled values. 
The concept of conventional AGC is discussed in [l]-[3]. 
Controlling the frequency has always been a major subject in 
electrical power system operation and is becoming much more 
significant recently with increasing size, changing structure 
and complexity in interconnected power systems. In a 
competitive electricity market, there will be many market 
players, such as generating companies (Gencos), distribution 
companies (Discos), transmission company (Transco), and 
system operator (SO). For stable and secure operation of a 
power system, the SO has to provide a number of ancillary 
services. One of the ancillary services is the "frequency 
regulation" based on the concept of the load frequency 
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control. A detailed discussion on load frequency control 
issues in power system operation after deregulation is 
reported in [4] -[8]. The AGC in a restructured electricity market 
should be designed to consider different types of possible 
transactions such as PoolCo-based transactions, bilateral 
transactions, and a combination of these two [7]. In this new 
paradigm, a Disco can contract individually with a GENCO 
for power and these transactions are done under the 
supervision of the System Operator (SO). Hence, situations 
may arise where DISCO make contract with only hydro 
generating units in competitive market place. In this paper 
we have analyzed the issue of AGC in restructured power 
system where contracts of DISOC exist only with GENCO 
based on hydro generating units. 

II. LINEARIZED MODEL OF AN INTERCONNECTED TWO- AREA MULTI 
UNITS HYDRO-HYDRO POWER SYSTEM 

Fig. 1 shows the schematic of an interconnected two-area 
power system. Each area has two GENCOs based on hydro 
generation andtwoDISCOs. LetGENCOl, GENC02, DISCOl, 
and DISC02 be in area I and GENC03, GENC04, DISC03, 
and DISC04 be in area II. Fig. 2 shows a complete block 
diagram of an interconnected two-area multi-units all-hydro 
power system for AGC study. The nominal system parameters 
are given in appendix [9]. The closed loop system in Fig. 2 is 
characterized in state space form as 

x = A cl x+B cl u 

where x is the state vector and u is the vector of power 
demands. State matrix ^i is constructed from Fig. 2 and given 

in appendix. Eigenvalue analysis of the state matrix j^i for a 
given operating conditions shows the positive real parts for 
some eigenvalue pairs (shown as shaded in Table 1 ) which 
indicate that the system itself is inherently unstable. Thus, 
load frequency stabilization can never be possible for such 
an all-hydro test system after the load perturbation. 
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Figure 1. 



Schematic of two-area system in restructured power 
system 
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Figure 2. Linearized model of an interconnected two-area restructured power system 
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Figura-4. Structure oFTCPS as frequency stabilizer 
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Figure 3. Area frequency response of system after 10% step load 
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Eigenvalue Analysis of Test System 
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Fig. 3 shows the transient response of the system after 
10% step load perturbation in area 1. It has been observed 
that under the occurrence of load changes, a system fre- 
quency is heavily perturbed from its normal operating point 
and system does not regain its stable state. To overcome this 
inherently unstable situation and to make the system stable 
against any load disturbance, the impact of Thyristor Con- 
trolled Phase Shifter (TCPS) [10] and Superconducting Mag- 
netic Energy Storage (SMES) [11] are analyzed. The struc- 
ture of TCPS and SMES as frequency stabilizers are shown 
in Figs. 4 and 5 respectively. TCPS is located in series with 
tie-line near area lwhereas SMES is placed at the terminal of 
area(s). 
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Figure 4. Structure of SMES as frequency stabilizer 
III. MATHEMATICAL PROBLEM FORMULATION 

The objectives of AGC are to reestablish primary 
frequency regulation, restore the frequency to its nominal 
value as quickly as possible and minimize the tie-line power 
flow oscillations between neighboring control areas. In order 

to satisfy the above requirements, gains (K n , K G ) of integral 
controller in AGC loop, parameters of TCPS (K TCPS ,T TCPS ) and 

the parameters of SMES (K SMES , T SMES ^.T^T, and T 4 ) ) are to 

be optimized. In the present work, an Integral Square Error 
(ISE) criterion is used to minimize the objective function 
defined as: 

Figure of Demerit(FDM)= V Im? +Af 2 2 +AP t - e 1 AT (1) 

where AT = a given time interval for taking samples, A.f = 

incremental change in frequency and AP tie = incremental 

change in tie-line power. The objective function is minimized 
with the help CRPSO based optimization techniques. 

IV. Craziness Based Particle Swarm Optimization 

The PSO was first introduced by Kennedy and Eberhart 
[12]. It is an evolutionary computational model, a stochastic 
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search technique based on swarm intelligence. 
Velocity updating equation: 



Jfc+1 

V: = V 



•■ +clxrlx(pBest ; -xM+c2xr2x(gBest-xf ) 



(2) 



Position updating equation: 

Jc f+i=xf+vf+ 1 (3) 

The following modifications in velocity help to enhance the 

global search ability of PSO algorithm as observed in CRPSO 

[B]. 

(i) Velocity updating as proposed in [13] maybe stated as in 

the following equation: 



r2xy*+(l-r2)xclxrlx(pBest -x ) 
+ (l-r2)xc2x(l-rl)x(gBest-x ) 



(4) 



Local and global searches are balanced by random number r2 
as stated in (5). Change in the direction in velocity may be 
modeled as given in the following equation: 

v- +l = r2xsign(r3)xvf +(l-r2)xclxrl 

xfpBest,. - xM+ (1 - r2)x c2 x (l - rl)x (gBest - xf ) (5) 

In (5), sign (r3) maybe defined as 



sign(r3) = 



-1 (r3<0.05) 
1 (r3>0.05) 



(ii) Inclusion of craziness: Diversity in the direction of birds 
flocking or fish schooling maybe handled in traditional PSO 
with a predefined craziness probability. The particles may be 
crazed in accordance with the following equation before 
updating its position. 



,/t+i 



,/;+! 



: v- • + Pr(r4) x sign{rA) x v f™"™ (6) 

where, Pr(r4) and sign(r4) are defined respectively as: 



Pr(r4) = 



1 (r4<P cm ) 
{r4>P aaz ) 



(7) 



sign(r4) 



(8) 



1 (M>0.5) 
-1 (r4<0.5) 

During the simulation of PSO, certain parameters require 
proper selection, as PSO is much sensitive to the selection of 
input parameters. The best chosen maximum population size 
= 50, maximum allowed iteration cycles = 100, best Pcraz = 0.2 
(chosen after several experiments), best values of cl and c2 
are 1.65. The choice of cl, c2 are very much vulnerable for 

PSO execution. The value of v . cmzmess lies between 0.25 and 

0.35. The novelty of CRPSO lies with the fact of its faster 
convergence to the true optimal solution as compared to its 
other counter parts. 

V Simulation Results And Discussion 

The concept of a "DISCO participation matrix" (DPM) is 
used to make the easier visualization of contracts between 
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GENCOs and DISCOs [7]-[8]. Each area is having twoGENCOs 
and two DISCOs. Let GENCOl, GENC02, DISCO 1 and 
DISC02 be in areal and GENC03, GENC04, DISC03 and 
DISC04 be in area2. Unlike in the traditional AGC system, a 
DISCO asks/demands a particular GENCO or GENCOs for 
load power. Thus, as a particular set of GENCOs are supposed 
to follow the load demanded by a DISCO, information signals 
must flow from a DISCO to a particular GENCO specifying 
corresponding demands. 

The demands are specified by cpf (elements of DPM) and 

the pu MWload of a DISCO. These signals which were absent 
in traditional AGC will carry information as to which GENCO 
has to follow a load demanded by which DISCO. 

apf l (i = 1,2,3,4) are the Area Control Error (ACE) 

participation factors of different GENCOS. 



DPM = 



cpfU cpf 12 cpf 13 cpfU 

cpfll cpfll cp/23 cpflA 

cpf 3 1 cpf 32 cpf 33 cpf 34 

cp/41 cp/42 cpf43 cp/44 



(9) 



For TCPS-SMES coordination, TCPS is considered in 
series with a tie-line near area 1 and SMES is placed at the 
terminal of area2 whereas for SMES-SMES coordination, 
SMES units are considered at the terminal of both areas. 

Case 1: PoolCo based transaction: 

An event is simulated in which a system shown in Fig. 2 
is subjected to step load perturbations as given in (10) [7] -[8] 
at t = 5 sec. 

The GENCOs in each area participate equally in AGC, indi- 
cated by Area Control Error (ACE) participation factors 
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(10) 



qfl = 0.5. q^j = 0.5. £^3 = 9.5. qf± = 0.5 



The 0.05 pu load is demanded by DISCO 1 and DISC02 
identically from their local GENCOs in areal only, i.e. GENCOl 
and GENC02, which is indicated by each non-zero 
participation factor of 0.5 in DPM. Similarly, DISC03 and 
DISC04 demand 0.05 pu load from their local GENCOs in area 
2, i.e. GENC03 and GENC04. As there is no contracts between 
inter-area GENCOs-DISCOs, the inter-area participation 
factors of DPM (shown shaded in (10)) are zero. This type of 
overall transaction is called as PoolCo transaction. The un- 

contracted load demand Af Luncotl and A/^ uncot2 of the areal 

and area2 respectively are zero. The transient response of 
the system for the transaction given by (10) is shown in Fig. 
6. It is clearly observed that the coordination of SMES-SMES 
and TCPS-SMES effectively suppress the area frequency and 
tie-line power oscillations and stabilize the system effectively. 
The optimized parameters for both coordinated system are 
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given in Table II. Power generation response shown in Fig. 6 
indicates that the initial surge of power is in the opposite 
direction to that desired due to non-minimum phase 



characteristic of hydro turbine [9]. Hence, issue of LFC with 
only hydro generating unit has become very critical and needs 

a 



— g 












SkES-SUE 

TPS-SHB 




d -a< 













_ an; 

™ o 
= -oik 



■o 2a ao tn 50 




■c so ao 




'o a a -o so 

~ r '< 3< : 





■ — ; : : 




P^~~ 


i i — -n 




L 


-n 


S'JBS-S'.'S 

gn nrn 








■ ' 






ri-p^ — i — : — 

jJ-flS ■ - -WT r r ! 

b ■ W ! ! ! ! 

m. ' ~ ~ H ~ I I I ----- 



: ■: :: :; -: 







o — -- */-r^2* i l 

-Qa - - Ht"-- j , ! r 

-■ . _ J - T , 1 K 

If I I ' ' 



—as 

± c 

t-GS 



L . § r 

ii 11 

1 i ' ! . 



3C 3C 



1 


1 1 1 


1 J 












Figure 5. Case 1 



Figure 6. Case 2 



Figure 7. Case 3 



special coordinated control of TCPS-SMES or SMES-SMES. 
Since there are no contracts of power between a GENCO in 
one area to a DISCO in another area, the scheduled steady 
state power flow over the tie line is zero. 

Case 2: Bilateral Contracts: 

The nonzero participation factor of DPM of (11) repre- 
sents both PoolCo and bilateral types of transactions unlike 
that indicated by DPM of (10). The apf, of GENCOl is the 
highest; the apf 2 of GENC02 is the lowest which indicates 
that the participations of GENCOl and GENC02 in zeroing 
ACE are the highest and the lowest respectively. The total 
generation required of individual GENCOs can be calculated 
as [8]: 
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The mutual scheduled tie-line power flows among the areas 
can be represented as [8]: 
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(11) 



(12) 



SchsiAdsd _ P t jj j_7 = 

The dynamic responses for combined PoolCo and bilateral 
transactions are shown in Fig. 7. In this case also, the coordi- 
nation of SMES-SMES and TCPS-SMES 

Table U 

Optimized Parameters obtained by CRPSO for TCPS-SMEs 

and SMES-SMES Coordination 
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works satisfactorily. As the inter-area bilateral contacts exist 
between the DISCOs and GENCOs, scheduled tie-line power 
are observed on the tie-line in steady state also. 

Case 3: Contract violation 

It may happen that a DISCO violates a contract by 
demanding more power than that specified in the contract. 
This excess power is not contracted out to any GENCO. This 
un-contracted power must be supplied by the GENCOs in 
the same area as the DISCO. It must be reflected as a local 
load of the area but not as the contract demand. Consider 
case 2 again with a modification that DISCO demands 0.05 pu 
MW of un-contracted excess power in (11) i.e. 

AP L,uncoti = 0.05pu . Fig. 8 shows the stabilized transient 
response in the event of contract violation with the 
coordinated device operations. 

The convergence characteristic of the objective function 
given in (1) with the help of CRPSO algorithm is shown in 
Fig. 9. 
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VI. Conclusion 

Eigenvalue analysis of an interconnected two area multiple 
unit hydro-hydro power system results in positive real parts 
for some of the eigenvalue pairs which confirms its 
dynamically unstable behavior. The optimized coordination 
of SMES-SMES and TCPS-SMES with the help of CRPSO 
effectively stabilize the system and significantly improves 
the transient response of area frequency and tie-line power 
exchanges under different contract variation in restructured 
electricity market. The performance of SMES-SMES 
coordination yields the least value of Figure of Demerit and 
shows lesser undershoot and overshoot in dynamic 
responses of frequency and tie-line power exchanges as 
compared to TCPS-SMES coordination. 
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